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Abstract 

We investigate the cooling rate of a gas of inelasticaUy interacting particles. 

When we assume velocity dependent coefhcients of restitution the material 

cools down slower than with constant restitution. This behavior might have 

large influence to clustering and structure formation processes. 
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The behavior of granular gases has been of large scientific interest in recent time. Gold- 
hirsch and Zanetti [|T| and McNamara and Young have shown that a homogeneous gran- 
ular gas is unstable. After some time one observes dense regions (clusters) and voids. To 
evaluate the loss of mechanical energy due to collisions one introduces the coefficient of 
(normal) restitution 

9' = eg , (1) 



where g = \g\ and g' = g , describing the loss of relative normal velocity g' of a pair of 
colliding particles after the collision with respect to the impact velocity g 

It can be shown that even for three particles for certain region of the coefficient of 
restitution there exist initial conditions which lead to a behavior which is called "inelastic 
collapse" . This means that the particles accomplish an infinite number of collisions in finite 
time 1^. The conditions under which one can observe inelastic collapse have been studied in 
one dimensional systems as well as in higher dimension |Q. Recently it was shown numer- 
ically that the probability for a collapse rises significantly when the particles have rotational 
degree of freedom [^. In this case the collapse is possible for much larger coefficients of 
restitution than for non rotating particles. Other interesting related results concern bounc- 
ing ball experiments on vibrating tables where complicated dynamical behavior is observed 
(e.g. 0). Recently complicated, under certain circumstances irregular motion of a bouncing 
cantilever of an atomic force microscope when excited by a transducer was investigated [0 . 

In the investigations the approximation of constant coefficient of restitution was 
assumed. Solving viscoelastic equations for spheres currently it was shown that the coef- 
ficient of normal restitution e is not a constant but a function of the impact velocity e {g) 
itself For the "compression" ^ = _Ri + — |^i — ^2! of particles with radii Ri and R2 
at positions rl and r2 one finds 

i+p{e" + \A4ii) = Q (2) 



2 yVr^ 
3 m^ff (1 - z/2) 



P = TT-iff?^ — (3) 
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Y is the Young modulus, v the Poisson ratio and 

77111712 . 

meff = (4aj 

mi + m2 

are the effective radius and mass of the grains, respectively. A is a material constant de- 
pending on the Young modulus, the viscous constants and the Poisson ratio of the material. 
Equation (H) was derived under the precondition that the colliding spheres have impact ve- 
locity much less than speed of sound in the particle material. For details see P] . The initial 
conditions for solving are 

^(0) = (5a) 
m = 9 ■ (5b) 

The coefficient of restitution e of at time t = colliding spherical grains can be found from 
this equation relating the relative normal velocities = ^(0) at time of impact and at time 
tc, when the particles separate after the collision, i.e. tc is the collision time: 

e = -i (tc) li (0) . (6) 

The (numerical) integration of equation (|^) yields the coefficient of restitution as a function 
of the impact velocity (see fig. 1 in [Q) which is in good agreement with experimental 
data |TD[. Constant coefficient of restitution, however, does not agree with experimental 



experience |TT|]. Other theoretical work on this topic can be found e.g. in P^JTB[ |. 

Consider a gas of granular particles at a given initial granular Temperature Tq. Then 
the question arises how the temperature decreases with time due to inelastic collisions. This 



problem has been investigated earlier [ll4| , p!5| for the case of constant coefficient of restitution 



and the result is fs. also \lt 



T{t)=To{l+t/T)-' . (7) 

The time scale r is a material constant. The temperature decay (|^) is the origin of the 
cluster instabilities which have been investigated recently [|l],0|. 



The aim of the present paper is to derive an exphcite analytic expression for the coefficient 
of normal restitution e{g) as a function of the impact velocity g. A direct consequence of 
this result will be a refined expression for the temperature decay of a granular gas. 

The duration of collision for the undamped problem {A = 0) is given by [0: 



90 

t = ^- (8) 

We want to point out here that G° is a constant pure number, not depending on any material 
properties. Hence, depends only on the material constant p and on the initial velocity g. 
We use equation (||) to define a rescaled dimensionless time 0: 



Q = p^o gZ t (9) 



Using the abbreviations 



V = p^g (10a) 
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a = \A (10b) 



and a new set of variables 



Q = p!' g5 t = t (11a) 

x(e) = p'm (lib) 



we rewrite (Q) in the form 



X + av -'X^fx^v 5a;2 = (12) 



with X = -^x. We see that 



dx. . 1 dC . g 1 dx , . 

Hence the initial conditions in our new variables x and read 

x(0) = (14a) 

^{0) = m = v'^ (14b) 



Both equations of motion, (0) and (0) become special at x = or ^ = 0, respectively, i.e. 
all derivatives of third order and higher diverge. This will be shown for the case of x: 

d d / _i . — _2 3\ 

— — X = — 77^ 0;^^ ^x\/x + v 5a;2 
did rfB ^ ^ 

= av~'^Axy/x + -^^^j=\ -'^v^ix^ (15) 
y 2y/x ) 2 

Hence 

d^ 

lim — ttttX = ±oo . (16) 
x^o dQ^ ^ ' 

and so are the higher derivatives. Because of this singularity we must not expand x in 
powers of B. Because of the initial conditions x(B) has the form 

a;(B) =t;tB (l + r/(B)) (17) 
^(0) = , (18) 

which defines the function ^^(B). Using transformation (|l^ we find 



B17 + 2?) + av'-Q'ifi^l + 7]+ {av^\/Q + B2 j (1 + 7^)2 = . (19) 

In (|TD|) occur terms B'''^ and B^'^, therefore we expand rj in powers of -\/B 

00 

77 = ^ afc B^ (20) 

fc=0 

The first coefficient Oq vanishes because of the initial condition for x. When we require 

^=^ + a. + ... (21) 

to be finite at B = the second coefficient oi must vanish as well. With Taylor expansion 
of a/1 + ?7 and (1 + 77)2 for small rj we arrive at 

4 1^3 4 ^5 3 1 9 2^0 . . 

7] = -—av^e^ - — B2 + —av^oO^ + —a^v^e^ . . . (22) 
15 35 70 15 

and therefore 
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1 4 5 4 4 7 1 2 6 4 

x = f 5 B av c)2 ti5t)2 -I a v^'a 

15 35 15 



Rearranging the full series ( pSf ) one finds 

a; = t;! fe - —95 + — 6^ + . . . 
V 35 175 

= wtxo(6) + arai(e) + a^vix2{Q) . . . (24) 

f txo is the solution of the undamped (elastic) collision (s. dashed line in fig. 0). The full line 
in fig. |1| shows the damped motion according to eq. (pB]). The direct numerical integration 
of eq. (|T^) collapses with the full line. 

For X (|0c) where is the duration of the undamped collision one finds using (p^ : 



" (f ) = (f ) + """^ (f ) + "'"^"^ (?)+••• 

= t;^5o + avBi + a^yiSa + . . . , (25) 

which we do not need now but later on. 

Note that the coefficients Bk are constants, i.e. they do not depend on v nor on material 
constants. 

Equations (|^) and (12), respectively, hold for the entire collision. The collision starts 



with V and ends with v . For practical purposes we now define the term inverse collision. 
The inverse collision is a collision which starts at time 6c with relative velocity v and ends 
at time with relative velocity f , i.e. time runs in inverse direction during the inverse 
collision. The equation of motion for x*™, i.e. for a collision in inverse time follows from 
(|l^). Since the inverse collision starts with v we have to replace v hy v . Because of the 



time reversal we have to change the sign of time derivatives of odd orders, i.e. x —x 
The equation of motion for the inverse collision reads 



tnv 



6 



X 



^™ - a (v yh'^'^V^ + (w')-t (x'"") ^ = (26) 
A motion due to eq. (|26|) in normal time would be an accelerated one. However, we shall 



mention here that eqs. (|T^) and (|26|) describe strictly the same physical motion. The 
solution X*™ of the inverse problem can be derived from the solution of the direct problem 
replacing a — > —a and v v' . 

x^™(e') = {v')lxo{e') - av'xiie') + a\v')lx2{e') . . . (27) 

Now we determine the collision time Gc and the final velocity. One direct method to calculate 
©c would be to determine the solution of x{Q) = using Taylor expansion of x in the region 
close to e°. It can be seen easily that this method fails since all derivatives of with 
n > 3 diverge for O = G°. Therefore Be has to be calculated by an indirect methode. 
The problem will be subdivided into two parts (s. fig. |^): 

a) the motion of the particles x from B = to time B^ when x approaches its maximum 
and where x changes its sign, and 

b) from Bm to Be. 

In case of undamped motion where a = we have B^ = Be/2. In part b) we do not 
consider the collision itself but the inverse problem in the interval (B = 0, B^), with B^ 
being the time where x^"^"" approaches its maximum. The continuity of both parts means 
x(BJ =x-MBJ. 

For finite damping a 7^ we write B^ = Be/2 + 6 and B^ = (B^!)' /2 + 6' and remind 
that BO = (B")' . To get an expression for 6 we expand 

= (x«(f)+«o(f)+^|3-o(f )+...) 



and using xq = {vt xq is the solution of the undamped problem) 

1X1(90/2) 2N 



The expression ( |30D has to be inserted into the Taylor expansion of x (0^/2 + S): 



X 



(e»/...)..*(.(f).a„(f)4.„(f )....) 
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, /©n a^t;! X? (60/2) 



f 5 Xo -r^ + avxi 



2 J J 2 xo (60/2) 



+ a\lx2 ( ^ ) + ^ ("^ 



Hence 

4 .6/ /60\ Ix?r60/2U^ 



X (6„) = .1x0 (60/2) + avx, (6O/2) + ahl (x, - ^ 

Replacing again v —>■ v' and a — *• —a yields 

r, I ,aii (0c/2) „ / 2\ 

x-^ (6:„) = (.')^a:o (60/2) - av'x, (6O/2) + a^(.')t (x, - ^11111^ 

As explained above both solutions ( p3D and ( p5| ) have to be equal. With 

'60\ lx?(60/2) 



+ 



= X2 



2; 2xo(©°/2) 



we write 



.1x0 (^) + a.xi (^) + a^.l/? = (.')^xo (^) - a.'xi (^) + a\v')%^ 
We expand f ' in a 

f ' = f + awi + a^f 2 + . . . , 

and find 
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ysxo ( + avxi ) + a^f 5/? 



with 5v = avi + a^f2 + . . .. Writing (1 + in powers of a and comparing coefficients 
yields finally 

= f (l - at;^Ci + a^t;^C2 + • • •) , (40) 



with 



5-1 (f) 



Xi 1 


(f) 


Xo 1 


(f) 



C2 = ^{^r^\ . (41b) 

Since does depend neither on any material properties nor on the impact velocity g or 
V, respectively, Ci and C2 are pure numerical constants. Evaluating Ci and C2 in (^T]) 
numerically yields = 1.15344 and C2 = 0.79826. 
For coefficient of normal restitution one gets 



v' 



e = — = I- av^Ci + a'^v5C2 + . . . (42a) 



V 



1 - CiAp^ + C2A^ p5 g-> + . . . (42b) 



with g being the impact velocity (fig For the duration of collision we find with (^ 
and (M^ 

= QcP'h (1 + ^C.aplg'^^ + 0{a') . (43) 
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9c = f 5tr 



= e°(^l + ^Ciat;^)+0(a2) (44) 

To check the theoretical result (eqs. ^1]) we integrated numerically eq. (|T2|) and received 
the curves e(f ) and 6c(f )• Then we fitted Ci and C2 to these data using (^) and (|^). For 
instance for a = 0.05 we found = 1.15356 and C^'''^ = 0.80439 from the curve e{v) 

(s. eq. pp. The fit of Ci to Qc{v) (s. eq. pi) gives Cf""* = 1.15342. For other values of a 



we found very similar numbers. Hence, the numerical results agree with theory. 

When we use the velocity dependent coefficient of restitution in the collision term of 
Boltzmann equation 

T ~ 1 1 dvidv2 (1 - e^) \v, -V2ff (v,) f {V2) (45) 
we get the cooling rate for dissipative gas 

5 

T ~ To/ (1 + :;7) ' (46) 

Our final results eq. (^2|) shows that for viscoelastic colliding smooth bodies the coefficient 
of normal restitution is a decreasing function with rising impact velocity: 1 — e ~ (7^. A 
direct consequence is the cooling rate of a granular gas (eq. (^B[)): a granular gas consisting 
of viscoelastic particles cools down significantly slower than a gas of particles which collide 
with constant coefficient of restitution (s. eq. (|^)). Due to our understanding it is neither 
self-evident whether the clustering observed in granular gases of the latter type, and the 
extreme case of this effect, the inelastic collapse, will change their overall behavior nor 
whether they exist at all. These questions should be reconsidered in detail for velocity 
dependent restitution. 

The authors are grateful to N. Brilliantov, S. Esipov, H. Herrmann, F. Spahn and 
W. Young for helpful discussions, and J.-M. Hertzsch for providing relevant literature. 
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FIG. 1. The dynamics of the coUision. The dashed hne shows the (strictly symmetric) solution 
of the undamped collision. For the case of the damped motion (full line) the maximum penetration 
depth is achieved earlier whereas the duration of the collision is longer (Og > 0^)- 
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FIG. 2. Sketch of the calculation. The first part (a) 6 G (0, 0^) is calculated directly, for the 
other part (6) we define the inverse collision where the particles start with velocity v' and velocity 
approaches zero at B = Qm- Both curves have to fit together smoothly. 
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FIG. 3. The coefficient of restitution over impact velocity due to eq. (p^. As expected for 
small relative velocity the particles collide almost elastically. The result of numerical integration 
of (P) coincides with the curve. Both curves cannot be distinguished in the plot. 
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